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Abstract 



We construct cubature formulas on spheres supported by homothetic im- 
ages of shells in some Euclidian lattices. Our analysis of these cubature for- 
mulas uses results from the theory of modular forms. Examples are worked 
out on S"- 1 for n = 4, 8, 12, 14, 16, 20, 23, and 24, and the sizes of the 
cubature formulas we obtain are compared with the lower bounds given by 
Linear Programming. 

1 Introduction 

For a dimension n 2 and a positive number r, let §™ _1 denote the sphere of 
equation (x \ x) = r in the Euclidean space R n given with the canonical scalar 
product (• | •), and let a denote the rotation-invariant measure on such a sphere, 
normalised by cr(§™ -1 ) = 1; we write S n_1 for §™ _1 . For an integer t ^ 0, a cubature 
formula of strength t on is a pair (A, W), where AT is a finite subset of S™ -1 

and where W : — * K>o is a positive- valued function such that 



for every polynomial function / : K" — » C of degree at most t. The set A is the 
support of the cubature formula {X, W), the points in X are its nodes, the values 
W(x) are its weights, and the cardinality \X\ of A is its size. A spherical t-design 
on is a finite subset X of S™ -1 which, together with the constant weight (by 
necessity of value \X\~ ), is a cubature formula of strength t. 

There are existence results according to which, for any dimension n and strength 
t, there exist cubature formulas, and indeed spherical designs, of strength t on S" _1 . 
For cubature formulas, there are elementary arguments using general methods of 
convexity, which give some upper bounds on sizes (see Section the existence 
of spherical designs is a particular case of results of ScyZas84 , which do not 
provide meaningful size bounds. There are no known straightforward constructions, 
except for n = 2, in which case, for any t ^ 0, the (t + l)th roots of unity in C 
provide a spherical f-design on the circle; and in larger dimensions for low values 
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of t (antipodal pairs, regular simplices, and regular hyperoctahedra are respectively 
spherical 1-, 2-, and 3-designs, of the lowest possible sizes). 

It is therefore a natural question to ask for more explicit constructions. Some 
have been given in terms of orbits of finite groups on the sphere (see |HarPac04| ) . 
and by other various methods (see |Bajn91| and later papers by Bajnok, |Kupe-a| , 
|Kupe-b| , and |HarlSlo96| for § 2 ). 

The object of the present paper is to describe constructions involving shells (or 
layers) of integral lattices in R ra . They are elaborations of constructions which can 
also be found in [Mart VOT] . |Pach-a| . and |HarPac-a| . 

Let E be a discrete subset of M™. For any m > 0, denote by E m the shell 
ED . If mi, . • • , m r are pairwise distinct positive numbers such that each shell 
E mj is non-empty the union of the sets ^1— E mj is a finite subset (possibly with 

multiplicities) of the unit sphere in M". Spherical designs and cubature formulas 
which appear below have supports of this kind, where E is either an integral lattice 
in K n or a union of lattices; observe that they are all antipodal, namely that X = —X 
and W(x) = W(—x) for all x G X. 

Results on S 3 Let us describe some examples provided by this method when 
n = 4, and compare them to known lower bounds for sizes and to other known 
examples. For n J? 4, the root lattice of type D„ is the integral lattice 

D n = {xeZ n \xi + h x n = (mod 2)}. 

(We denote root systems as D„ by bold letters and lattices as D n by thin letters.) 
It is important for our computations that Z3 4 is a 2-modular lattice (the definition 
is recalled in Section 21 below). 

For n 5, n ^ 7, it is known that each shell {D n )2m is a spherical 3-design 
which is in general not of strength 4; for n = 4, each shell (Di)2m is a spherical 
5-design which is in general not of strength 6 (it has been checked in |Pach-a| that 
it is never of strength 6 for 2m ^ 1200). The sizes of the shells of the root lattice 
of type D4 are given by the coefficients of the theta series 

e Di = J2 \( D ihm\ q 2m = n(2m)q 2m 

= 1 + 24 q 2 + 24 q 4 + 96 q 6 + 24 q 8 + 144 q 10 + 96 q 12 + 192 q 14 + ■ ■ ■ 

where ri{2m) denotes the number of ways of writing 2m as a sum of four squares. 

The dual D\ = {y e W 1 \ (y | D4) C Z} is similar to D4 and the renormalized 
dual D' 4 — \[2D\ is isometric to D4; moreover the intersection D4 D D' 4 is reduced 
to {0}, since the non-zero coordinates of the vertices in D 4 are never in Z. For our 
analysis, we use some invariant theory of the finite group Aut(i?4 U D 4 ). 

Theorem D 4 . Set E = L> 4 u D 4 

(i) For any even integer 2m ^ 2, the shell E 2m is a spherical 7-design; in partic- 
ular, E2 is a spherical 7-design of size 48 in S 3 - 

(ii) The two shells £2, E§ provide the support X — ^-^2 U ^7g^6 of a cubature 
formula of strength 11 and size 240. 

(Hi) Similarly, -±E 2 U ±£ 6 U -^E 10 is the support of a cubature formula of 
strength 15 and size 528. 

Claim (i) for E 2 appears already in Section 6 of |GoeSei8l] . where the authors 
observe that E 2 is an orbit of the Weyl group of type F4. Claim (i) of Theorem E s 
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below and the analogous result involving the two shortest shells of a Leech lattice 
appear also in |GoeS ei81 . 

For comparison, here are some known results concerning the 3-sphere. 

Dickson reports on a formula due to Liouvillc (1859) which essentially shows 
that the system of roots (.04)2 of type D4 is a spherical 5-design, a formula due 
to Kempner (1912) which essentially show that (D4 U ^4)2 is a spherical design 
of strength 7, and a formula due to Schur (1909) related to a cubature formula of 
strength 11 and size 144. See |Dickl9l pp. 717-724]. 

There is a formula of strengh 7 and size 46 described in HadSlo94 . 

There is an exceptional regular polytope in dimension four known as the 600- 
ccll. It has 120 vertices which constitute a spherical 11-design; it is an orbit of a 
Coxeter group of type H4. Any orbit of this group is a spherical design of strength 
11, but there is a special orbit of size 1440 which is of strength 19; see |GoeSei8Tl 
Section 5]. Another cubature formula of the same strength, 19, and of size 720 
appears in |Sali75j : its support is the union of the 120 vertices and of the 600 
centers of faces of the 600-cell. 

The best lower bounds known to us for the size of cubature formulas in S 3 are 
listed in Tables and of Section □ 

Results on § 7 Let i~g be a Korkine-Zolotareff lattice, namely an even unimodular 
lattice of dimension 8 (uniquely defined up to isometry by these properties). Such 
a lattice is generated by a root system of type E 8 . The sizes of its shells are given 
by the coefficients of the theta series 

©A = £ K^WI 3 2m = 240 ]T a 3 (m)q 2m 

m^O m^0 

= 1 + 240 q 2 + 2160 q 4 + 6720 q 6 + 17520 q s + ■ ■ ■ 

where cr 3 (m) = £ d | m d 3 . 

It is known that each shell (Is) 2m is a spherical 7-design which is in general 
not of strength 8. Indeed, the shell (r 8 ) 2m would be a design of strength 8 if and 
only if the Ramanujan function t would vanish at m; it is a conjecture of Lchmcr 
|Lehm47j that r(m) ^ for all m ^ 1, and this has been checked for m ^ 10 15 
|Serr85| : more on this in |Pach-aj and in Item 1.16 of jHarPac-a] . 

Theorem E s . Set E = f 8 . 

(i) The two shells E2, E4 provide the support X = E2 U \E± of a cubature 
formula of strength 11 and size 2400. 

(ii) Similarly, -^E 2 U 4j-2'6 U E% is the support of a cubature formula of 
strength 13 and size 24240. 

(Hi) Similarly, -^E 2 U \E± U -^E 6 U E$ is the support of a cubature formula 
of strength 15 and size 26400. 

The largest known lower bounds for the sizes of cubature formulas on § 7 are 
given in Table |H1 of Section 

Dimensions 12, 14, 16, 20, 23, and 24 There are similar constructions on S 23 
in terms of a Leech lattice, an even unimodular lattice A of dimension 24 with empty 
shell A2 (uniquely defined up to isometry be these properties). We also describe 
constructions with shells of 

• the 3-modular Coxeter- Todd lattice K12, 
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• the 3- modular Quebbemann lattice Q14, 

• the 2-modular Barnes- Wall lattice BWiq, 

• the 2-modular Ncbc lattices iV 2 o (three of them), 

• the 3-modular Nebe lattice N 2 4 7 and 

• the shorter Leech lattice O23 (which is up to isometry the unique unimodular 



integral lattice A in R 23 with A x = A 2 = 0). 

Constructions with one shell of these lattices are already described in |BacVenOT] . 
Our results appear in Tables |5] to [21 of Section 

2 Bounds on the size of cubature formulas 

For fixed values of n and t, we are interested in constructing cubature formulas on 
of strength t of small size; it is therefore useful to estimate the minimal size 
that such a cubature formula can have. We recall in this paragraph some general 
results on this subject. 

Let b(n, t) be the minimal size of a spherical cubature formula on S™" 1 of 
strength t. 

Linear Programming bound [D eGoSe77L Theorem 5.10], see also [Yudi97 

Let Q^ k ', k S Zj,o> be the orthogonal polynomials for the scalar product 

(f\g) ■■= r f{t)g{t){i-t 2 )^-^dt 



normalized to (Q^* 1 | Q^ fc ') = Qw(l). For a continuous function F : [—1, 1] 
we define the numbers Ft G C by 

_ (F I Q( fc >) 
Pk Q(*)(l) ' 

so that F = £ fc>0 F k QW. Set 



M(n,t) := \F : [-1, 1] — ► R continuous 



F(u) ^ for u E [-1,1], 
F^O, and F fc < for k > tj' 

then b(n,t) ^(1)/^ for every F £ M(n,t). In other words 

F(l) 

b L p(n,t) := sup — — 

FeM(n.t) 

is a lower bound for b(n, i). Moreover, it is a consequence of |NikNik02l Therorems I 
and 2.7] that there exists a polynomial function F S M(n,t) for which bLp(ra,t) = 
F(l)/F . 

When i is odd, we show: 

1. Proposition. Let t = 2s+ 1, and Jet 

M{n,t) := {F e M(n,t) | F(-tt) = F(u)}. 

Then 

, / x 2F (!) 
b LP (n,t)= sup — - — . 

FeM{n,t) 
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Proof. Let F G M(n,t); set F(u) = + F(-it)) /2; we have F G M{n,t). Set 

F*(u) = (1 + u) One checks that 

F* G M(n,t). 

Yet 2F(1) = F(l) + F(-l) ^ F(l) and F = F = F *, therefore 

F(l) 2F(1) _ F*(l) 
~F^^ F ""^T' 

The result follows. □ 



In general, the exact value of bL,p(n, t) is not known. The following bounds are 
obtained by chosing a particular F in M(n,t): 

Delsarte (or Fisher-type) bound [DeGoSe77, Theorems 5.11 and 5.12] 

Choose F(u) := (ELp Q W K>) 2 if * = 2«, andF(u) := (l+u) (E?io Q (i_!y) («)) 2 
if t = 2s + 1. We obtain the bounds 

b D (n,2s+ 1) = 2 

Moreover, it is known that, if n ^ 3, the equality b(n,t) = br>(n,t) is possible 
only for some values of (n,t): see BanDani79]. |Ean Dam80 . and BaM uVe02] . 
There also exists a criterion to decide when bn{n,t) — b^p(n,t): see Ni kNik02l 
Theorem I]. 

Yudin bound [Yudi97 Another choice of F gives the bound 

£ Jl - u 2 )(»- 3 )/ 2 du 
by(n,t) = — = , 

J^(l _ U 2)(n~3)/2 du 

where 7 is the largest root of the polynomial (Q'* +1 ')'. 

Special cases For some values of (n, t), it is possible to compute the exact value 
of b(n,t). For example, it is shown in |AndrOO| that b(4, 11) = b LP (4, 11) = 120. 
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Numerical estimate of the Linear Programming bound For fixed values of 
(n, t), t = 2s + 1, the following procedure gives an estimate of the Linear Program- 
ming bound bLp(n,i). (There is a similar procedure when t is even.) 

• Choose a degree s, and a finite subset A C [0, 1] of well-distributed points, 
for example A = {i/N \ i = 0, 1, . . . , N} with N large. 

• By linear programming, find, among all polynomial F satisfying 

d 

F = 1 + FkQ (2k) , F(u) ^ for u G A, F k < for k > t, 

i=l 

the polynomial minimizing F(l); denote it by G. 



■5 



• Set e := — inf x g[o,i] G(x) ^ 0. The polynomial G := G+e is in the set M(n, t), 
and we have the estimate 

2 6(1) 2(G(l) + e) , , , 

-J-> . -L_J <bLp(nii) . 

Note that we have also b^p^t) ^ 2(7(1), where 

b$(n,t):= sup 

FeM(n.t) ^0 

Yet, if ci is sufficiently large and A sufficiently dense, 2G(1)/Gq is a good approxi- 
mation of t>Lp(?i,i)- In practice, we first apply the procedure with a relatively large 
d and a relatively small set A; we observe that Gk = for k > do- Then we apply 
the procedure with d = do and with a larger set A. 

The procedure just described imitates the one used in |ConSlo99l Chap. 13] to 
compute bounds for kissing numbers. 

We end this section with another kind of bound. Set 

fn + t-VS fn + t-i 

B(n,i) := 



n — 1 7 V n — 1 



Note that B(n,i) = bo(?i, 2i). This number is the dimension of the space of the 
restrictions to the sphere of the polynomial functions of degree at most t. 

2. Proposition. For every n and t, there exist cubature formulas with at most 
max{l,B(n,i) — 1} nodes. 

3. Proposition. Let (X, W) be a cubature formula on S" -1 of strength t. Then 
there exist a subset X' C X and a weight function W' : X' — > R>o such that 
\X'\ ^ B(n, t) and (X 1 , W') is a cubature formula of degree t. 

See |HarPac-al Proposition 2.6]. The second proposition is a consequence of 
the proof of the first one and provides an algorithm for computing (X' ', W') from 
(X,W). 

Note however that these propositions do not hold for spherical designs (case of 
constant weights). 



3 Harmonic polynomials 

It is sufficient to check condition (CF) in the definition of cubature formulas for 
harmonic polynomials only. Recall that a smooth function / on K™ is harmonic 
if A/ = 0, where A = J^^ii^/dxi) 2 is the usual Laplacian. We denote by 
J{( k ) (R n ) the space of harmonic polynomial functions M™ — > C that are homogenous 
of degree k. It is a classical fact that the space of restrictions to S™ -1 of homogenous 
polynomials of degree at most t coincides with the space of restrictions of the 
direct sum ® fe=0 7i^(W n ). This implies the following well-known criterion; see for 
example VenMarOl, Theorem 3.2] (details are written there for spherical designs, 
but the proof carries over to cubature formulas). 

4. Proposition. Let X be a finite subset of some sphere S^," 1 in E n , let W : X — > 
R>o be a weight function, and let t ^ be an integer. Then (X, W) is a cubature 
formula of strenght t on S"^ 1 if and only ifj^xex W{x) — 1 and 

W(x) P(x) = for all P € H (k) (M. n ), 1 < k < t. 
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Let 0(n) be the orthogonal group of the Euclidean space 1"; this group acts 
by isometries on the sphere S™ -1 . Moreover, O(n) acts naturally on Tt^ (R™) by 
.9 ■ / = /°.9 _1 - 

Whenever a group G acts on a space V, we denote by V the subspace of all 
elements of V fixed by G. 

Proposition^] above can be refined in the case the cubature formula is invariant 
by a finite subgroup of G. Namely, we have the following result, known as Sobolev's 
theorem ( |Sobo62| . see also |Sobo96l Chap. 2, § 2, Theorem 2.3]): 

5. Proposition. Let X be a finite subset of S" -1 and W : X — > R >0 be a weight 
function. Let G be a finite subgroup of 0(n) leaving (X, W) invariant, that is 

jiEl andW{gx) = W(x) Va; <E X,Vg G G. 

Then (X, W) is a cubature formula of strength t if and only if 

W(x) P(x) = VPeW (fe) (M n ) G ,l < 

xex 

Proof. For P £ H (k) (R n ), let P G := \G\~ 1 Y, 9 &g9 ' /• 0n the one hand > the 
G-invariance of (X, W) implies 

Y;W{x)p g {x) = y. w ^ p ( x ^ 

xex xex 
and, on the other hand, the G-invariance of the Lebesgue measure implies 

/ P G (x)d<r(x)= P{x)da(x). 
The result follows from Proposition 01 □ 



4 Modular forms 

This section is a reminder from [BaqEnOl | ■ Consider a lattice A in R n and a 
homogeneous harmonic polyomial P S TC^(M. n ) of degree k. The theta series of A 
with harmonic coefficient P is the formal power series 0a, p defined by 

x£A m>0 ^igA,„ ' 

In case P — 1, we write simply 0a- Set 

g := e 47rz , zeH, 

with H = {z e C | 3z > 0} the Poincare half-plane; the theta series Qa,p converges 
uniformly on every compact subset of H and defines consequently a holomorphic 
function on H. When the lattice A satisfies appropriate conditions, this function 
is a modular form of weight u), and we can use results on modular forms for the 
computation of the theta series. 

Let us first describe the class of lattices which plays the most important role 
below. A lattice A in 1" is even if (x | x) £ 2Z for all x £ A; such a lattice is 
contained in its dual A* = {y £ M™ | (A | y) C Z}, and there are integers £ ^ 1 such 
that \flA* is again even (£ is not unique since, if \/lk* is even, so is Vk 2 £A* for 
k ^ 1). For given integers n and £, denote by 

Cat n (£) 
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the class of even lattices A such that A' := \/£A* is even and det(A') = det(A); 
the latter condition implies that det(A) = l n l 2 . Clearly, if A € £at n {£), then 
A' e Cat n {£). A lattice A € Cat n {€) is l-modular if A' is equivalent to A. 

Now, we describe the modular forms associated to these lattices. For w a 
nonnegative integer and e e { + ,—}, we define 

as the space of holomorphic function / : H — > C that verify 

f(z+l) = f(z) 

'(-n)-(vc)"«-> 

for all z£i, and that are holomorphic at infinity, i.e., bounded on {z e H | Sz > 
yo} for y > 0. This means that / can be written as 



f(z) = a m q 2m , q = e 



ITVZ 



Let /(oo) := lim^^^oo f(z) = a . We define 

°Ml(l) = {feMl(e)\f(oo) = 0}. 

We give here the classical terminology: A modular form of weight ui > for a 
discrete group T C SL2(K) and a character \ '■ T — > K is a holomorphic function 
/ : H — > C that is holomorphic at infinity and that satisfy 



X((g5) ) ^ + 6 
(cz + c?) 



A parabolic form is a modular form / that is zero at infinity, namely such that 
/(oo) = 0. 

For £ a positive integer, let T*(f) be the subgroup of SL 2 (R) generated by 
r W = {(SS)eSL 2 (Z)| C =0mod^} and t t = ^f^j . 

For any integer s, let %s be the multiplicative character of T*(£) defined by 

x.((S3)) = (^) for(^)eroW, 

where ( — ) denotes the Kronecker symbol. Note that Xs depends only of the class 
of s modulo 4. Denote by Ai^(£) the space of modular forms of weight u> for the 
group r*(f) and the character Xw> an d by M~(£) the space of modular forms of 
weight u for the same group and the character x^+2 5 denote by °M^(£) the space 
of corresponding parabolic forms. 

Note that Equations (*) say that A4+(£), respectively A4~(£), is the space of 
modular forms of weight u for the subgroup generated by 

t =(i i) - *-(.°^ r 

and for the character \ u , respectively 
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6. Lemma. For £ £ {1, 2, 3}, the group T*(£) is generated by T and Si. in partic- 
ular, for these values of £, we have 

Ml{l) = MUZ) and °Mt(l)= Q Ml(£). 

Proof. Let I £ {1,2,3} and 7 £ { a c b d ) £ T*(t). We show, by induction on c 2 £ N, 
that 7 £ (T, St), where (T, Sg) denotes the group generated by T and Sg. 

If c = 0, then 7 = ( J ™ ) = T m for an integer m; thus 7 £ (T, S e ). 

If c ^ 0, let fc £ Z such that |d + kc\ < |c/2|. We have 



-(6 + *o)\/l a/V?\ 



and (-(d + kc)Vi~) 2 < c 2 (^/4) < c 2 . By induction, 7 T fc S' £ € (T,&), therefore 

y€(T,St). □ 

7. Proposition. Let A G £at n (£), where £ is a positive integer. Then 

e A + Q A , £M+ /2 {£), 

e A - e A / € °m- /2 (£). 

Let moreover P £ H (2h) (W l ), h^l. Then, if 2h = mod 4, 

Q^P + e A ^ P £°M+ /2+2h (£), 
®A,p-®A',pe°M- /2+2h (£), 

and, if 2h = 2 mod 4, 

e A , F + e A ,, P e °M- /2+2h (£), 

eA,p-eA',pe°M+ /2+2h (£). 

Partial proof. We give the proof only for £ £ {1,2,3}. The general case is more 
complicated: see |Ebel94l §3.1] together with jVenMarOll Chapter 2]. 1 

It is straightforward that the theta series involved are holomorphic in H and 
holomorphic at infinity. Now, by Lemma El it is sufficient to check that the theta 
series satisfy Equations (*). 

As A and A' are even, we have clearly Ql,p(z + 1) = 6l.p(z) for L = A or A'. 
Now, as a direct consequence of the Poisson Summation Formula — see for example 
|Ebel94l Prop. 3.1, p. 87] — we have, for any lattice L of rank n and any P £ 

e L ., P (z) = {fetL) l l 2 {-l) h {i/z) n ' 2+2h Q L , P {-llz). 
From this formula, we deduce, for L = A or A', 

®v,p{z) = (-l) h (i/V£z) n/2+2h e L , P (-l/£z). 
The result follows. □ 
Let 

M+(£) :=0M+W, 
°M + (£) := 0°M+(i), °M-(£) := 0°ACW- 

These are ./W + (£)-modules graded by the weight. For some values of £, the structure 
of these algebras is known: 

1 Note however that our examples of constructions of designs in Section [H] involve lattices in 
Cat n {£) for I g {1,2,3} only. 
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Remark. In general, the numerical subscript of modular forms given below indicates 
the double of its weight; for example A 2 4 S °.Mi 2 (l) is of weight 12, and 8e 8 € 
^4(1) is of weight 4. We have made an exception for the Eisenstein series Eha, for 
which we have retained the traditional notation. 

8. Theorem. Let £ be 1 or a prime number such that £+1 divides 24. Let k , k\ 
and fc 2 given by 

- 1, 

: 1 or 2 mod 4, and £ 7^ 1, 
: 3 mod 4, 




Then, we have 



where 



?2fc 



k 2 = k Q + k 1 + 2. 



.M + (£) = C[0 2fco ,A 2fel ], 
o M + (£)=A 2fel C[0 2 fc OJ A 2fel ], 
°^-(£) = $ 2fe2 C[0 2fco ,A 2fel ], 



9 Lo e M £,(£), with L € £at 2feo (£) ; 
fci 



A 2fel = (,7(z)^z)) M G X+W, 

and where Lo is given in TaWeQJ 

For a proof, see [rjacVenOTl Section 2] and |Queb95| Section 3] . Recall that 

Viz) = 9 1/12 ft (1 " <l 2m ) = 9 1/12 (1 - 9 2 - q 4 + 0(q 10 )), q = . 
m— 1 

For £ = 1, the modular forms of the theorem are A 2 4 = r)(z) 2i , whose Fourier 
coefficients are the Ramanujan numbers, 

©e 8 = L (4) , and $ 36 = A 24 E( 6 ), 

where E^ is the Eisenstein series of weight k. 

Tables to 01 give the theta series 2 fc o , A 2/ t 1 and <5> 2 fc 2 Ior au possible values 
of £. In the third column of Tabic ^ the lattice Lq is designated either by a root 
system that generates it (e.g., E 8 ), or by a quadratic form (e.g., ( 2 \))- In Tabled 
Pk is a suitable element of WW(K n ), the lattice L4 is the unique (up to isometry) 
4-dimensional lattice of minimum 4 and determinant ll 2 , and the lattices Lq 1 '' and 
L are the two lattices given in the column Lo of Table ^ We normalize $ 2 fc 2 by 
$2fc 2 = q 2 + 0(9 4 ). 



5 Construction of cubature formulas using shells of lattices 

Let G be a finite subgroup of O(n), and let X\ 1 X2,...,X r be nonempty finite 
subsets of § n_1 such that G Xj = Xj. Consider the set 

r 

X := (J Xj C S"" 1 . 
j'=i 
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i 


k 


Lo 


#2fc — ©Lo 


1 


4 


Eg 


1 + 240 9 2 + 2160g 4 + 6720g 6 + 0(q 8 ) 


2 


2 


D 4 


1 + 24g 2 + 24g 4 + 96q 6 + 0(q 8 ) 


3 


1 


A 2 «(?£) 


1 + 6q 2 + 6q 6 + 0(q 8 ) 


5 


2 


A 4 


1 + 20q 2 + 30q 4 + 60q 6 + 0(q 8 ) 


7 


1 


ill) 


l + 2q 2 + 4q 4 + 0(q 8 ) 


11 


1 


ill) 


1 + 2q 2 + Aq 6 + 0(q 8 ) 


23 


1 




1 + 2q 2 + 0(q 8 ) or 1 + 2q i + 2q 6 + 0(q 8 ) 



Table 1: Modular forms 02k o = ©l 01 weight ko 



I 


fci 


A 2kl = ( V (z) V (£z)) kl 


1 


12 


q 2 - 2Aq 4 + 252y 6 - U72q 8 + 0(q w ) 


2 


8 


q 2 - 8q 4 + I2q 6 + 64q 8 + 0(q w ) 


3 


6 


q 2 - 6q 4 + 9q e + 4q 8 + 0(q 10 ) 


5 


4 


q 2 - 4q 4 + 2q 6 + 8q 8 + 0(q 10 ) 


7 


3 


q 2 - iq A + 5q 6 - 7q 8 + 0(q 10 ) 


11 


2 


q 2 - 2q 4 -q 6 + 2q 8 + 0(q W ) 


23 


1 


q 2_ q i_ q6 + q^IO) 



Table 2: Modular forms A 2 fei = (v( z ) v{^ z )) 1 of weight ki 



I 


fc 2 


$2k 2 (definition) 


$2^:2 (expansion) 


1 


18 


©E s ,Pi4 


g 2 - 528g 4 - 4284g 6 + 147712g 8 + 0(q 10 ) 


2 


12 


©D^Pg - 0(Dl)',P 8 


q 2 - 88q 4 + 252y 6 + 64 g 8 + 0(q w ) 


3 


9 


@A1,P 6 + @(A1)',P 6 


q 2 - Uq 4 + A8q 6 + 68q 8 + 0{q 10 ) 


5 


8 


©Lo.Pe + ©L ,P 6 


q 2 - Uq 4 - A8q 6 + 68q 8 + 0{q 10 ) 


7 


7 


Qig.F, - (Lg)',P4 


q 2 - 10q 4 - Uq 6 + 68q 8 + 0(q 10 ) 


11 


5 


@L -LL 4 ,P2 + ®(Lo-LL 4 )',P2 


q 2 - 6q 4 - 3q 6 - Uq 8 + 0(q w ) 


23 


4 


e 4 1 )xL 2 ),P 2 + e (L< 1 )x4 2 ))',P 2 


q 2 - 2q 4 - 5q 6 - Aq 8 + 0(q 10 ) 



Table 3: Modular forms $2fe 2 01 weight k 2 
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The aim is to find numbers W\, W%, ■ • • , W r € R>o such that (X, W) is a cubature 
formula of high strength, where W is the weight function given by 

W(x)= E W i- 

j : xeXj 

By Proposition 01 (X, W) is of strength t if and only if the following conditions on 
W\ , . . . , W r are satisfied: 



3=1 

r 

E Wj P(Xj) = o, VP e w (fe) , l s$ fc s$ t, 

3=1 

where 

P{Xi) := £ W 

The following statement, together with Proposition \5\ is fundamental in our anal- 
ysis. 

9. Lemma. Let G be a finite subgroup of O(n), and let P C M" be a discrete set 
that is invariant by G. Let mi, TO2, ■ ■ ■ , w r be pairwise distinct positive numbers, 
let 

Xj = -^E m . C S"- 1 , ? = l,...,r, 



where E m = {x e £ | (1 | x) = to}. Let X = U£=i -^3 an d W\, . . . , W r > 0, and 
Jet k be a positive integer. 

Assume that there exist formal series 

Q i =^2a i (m)q m , i = l,...,N 

m>0 

such that, for every P S 7^*0 (R" ) G , the theta series &e,p is of the form 

Qe,p = ci{P) Q 1 + c 2 (P) e 3 + • • ■ + CJV (P) 9at 

for some Cj(P) G C. 

Then, for the condition 

r 

E WjPiXj) =0 VP e 7i (fc) (M") G 

3=1 

to LoJd, it suffices that 

«) E^^ = °> i = l JV. 

3=1 m 3 

Proo/. Assume that Equation ()}) holds, and let P <E W (fe) (R™) G . B Y hypothesis, 
there exist ci , • • • , cat € C such that 

E ^(^™)<z m = e^p = E c * e * = E ( E « a ^ m ) ) i m - 

m>0 i=l m>0 M=l ' 
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Since P is homogeneous of degree k, we have P((l/y/m)E m ) = mT k l 2 P(E„ 
Therefore, 



P(X i )=mT fc / 2 P(E m .) = ^, 



fc/2 



hence 



e 



A,P 



6 Applications 

If not otherwise specified, Lemma El is applied with the trivial group G = {id}. In 
this section, we do not give the details for all lattices we have considered, but we 
have made a selection which reflects most situations that may occur. All cubature 
formulas we have obtained are listed in Section below. 

6.1 The Leech lattice 

Let A be the Leech lattice, namely the unique (up to isometry) even unimodular 
lattice of dimension 24 and of minimum 4. Let P € H^ 2h ' (R™) . From Proposition 
and Theorem 03 we deduce (using A' = A): 

'e| g - 720A 24 ifP = l, 

if P e H (2/l '(M 24 ), 2/1 = 2,4,6,8,10,14, 

Cl (P)A 2 4 if Pe H( 12 >(R 24 ), 

.c 2 (P)e E8 A| 4 ifPGW(16)(R24 ); 

where C\ and c 2 are linear forms on H il2) (M. 24 ) and H (16) (R 24 ) respectively. The 
first coefficients of these theta series are 

9| 8 - 720 A 24 = 1 + 196 560<? 4 + 16 752 960 g 6 + 0(g 8 ), 
A 2 A = q 4 -A8q 6 + 0(q s ), 
9 Es A| 4 = + 192<? 6 + 0(g 8 ). 

Since 6 a ,p = for P G H (2h) (R 24 ) , 2 ^ 2h ^ 10, it follows from PropositionHthat 
every shell of the Leech lattice is a spherical 11-design. 

Now, we want to find cubature formulas of higher strengths using the construc- 
tion of Section [3] Consider the set 

X = Xi U X 2 , where X x = -^A 4 , A 2 = 4= A 6- 

V4 V6 

The two sets X\ and A 2 are disjoint. Indeed, if we had x G Xi fl A 2 , we would 
have y/Zx € A 4 and y/6x € Aa; but {y/Zx \ s/Ex) = 2y6{x \ x) = 2y/6, which is 
impossible since A is integral. 

We want to find out the numbers W\ and Wi such that (X, W) is a cubature 
of strength 15 (same notation as in Section |SJ). For this, we have to fulfill the 
conditions 



^2 W(x) P(x) = for P = H (12 \l 
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The first condition is equivalent to |A~i| W\ + \X 2 \ W2 = 1, that is 

196 560 W x + 16 752 960 W 2 = 1. 

By Lemma the second condition is equivalent to 

1 48 
— W-\ W-? = 

or W2 = (3/4) 5 Wi. The solution of this set of two linear equations is 

Wi ps 2.394 x 10^ 7 , W 2 ps 0.568 x 10~ 7 . 

So we get a cubature formula of size 196 560 + 16 752 960 = 16 949 520. 

We can obtain in the same way a cubature formula of strength 17 by using the 
shells of norm 4, 6, and 8, and a cubature formula of strength 19 by using the shells 
of norm 4, 6, 8, and 10 of the Leech lattice. 

6.2 The Korkine-Zolotareff lattice 

Let A be the Korkine-Zolotareff lattice, namely the unique (up to isometry) even 
unimodular lattice of dimension 8. It is the lattice generated by the root system of 
type E 8 . From Proposition [7| and Theorem HO we deduce (using A' = A): 

'e E8 ifP = l, 

if P eH( 2h \M. 8 ), 2/1 = 2,4,6,10, 

(->.v/> <( Cl (P)A 24 ifPeW^OR 8 ), 

C3 (P)6 E8 A 24 UPeH^H 

_c 2 (F)$ 36 \fPen^( 

So, by Proposition 0] any shell of the Korkine-Zolotareff lattice is a spherical 7- 
design, and we can obtain a cubature formula of stength 11 by combining the shells 
of norms 2 and 4. 

But if we try to obtain a cubature formula of strength 13 of nodes 

X = X 1 UX 2 UX 3 , where Xi = -=A 2 , X 2 = -J=A 4 , X 3 = —j=K, 

V2 v4 v6 

we obtain the weights 

W 1 ps -0.744 x 10~ 4 , W 2 ps 1.587 x 10" 4 , W 3 ps 1.005 x 10~ 4 , 

with W\ negative. A similar problem occurs if we try to use the shells of norm 2, 4 
and 8. Therefore, in order to get a true cubature formula with positive weights, we 
have to use 

X = X, U X 2 U X 3 , X l = -±=A 2 , X 2 = -J=A 6 , X 3 = i=A 8 . 

We obtain the weights 

Wi ps 0.792 x 10" 4 , W 2 ps 0.457 x 10~ 4 , W 3 ps 0.385 x 10~ 4 . 

Here, we have X\ C A3, since 2x £ A§ if x £ A 2 . Therefore, we obtain a cubature 
formula of strength 13 of size 

\X\ = \X 2 \ + |A 3 | = 6 720+ 17 520 = 24 240 

and of weights 

0.457 x 10~ 4 ifxeA 2 , 
0.385 x 10- 4 if x E X 3 \ Xi, 

W 3 ps 1.177 x 10- 4 if x G X\ . 




14 



6.3 The Barnes- Wall lattice of dimension 16 



Let A be the Barnes- Wall lattice of dimension 16. It is a 2-modular lattice, that 
is A' := \/2A is equivalent to A, and it is of minimum 4. From Proposition and 
Theorem [SJ we have: 



©a : p + ©A',p = 



2 9^ 4 -192Ai 6 ifP = l, 

if PeH^ h) {R 16 ), 2/i = 2,4,6,10, 

ci(P)A 2 6 if P 6H(«) (M 16 ), 

c 2 (P)A 2 6 6 2 54 ifPeW( 12 )(R 16 ), 

c 3 (P) $24 A 16 6 D4 if Pe7i( 14 )(M 16 ), 



@A,P — ©A' P = 



and therefore 



e 



A,P 





C4 (p) a 2 6 e D4 

C 5 (P) $ 24 Ai 6 

6^ -96A 16 


c 6 (P)A 2 6 
c 7 (P)A 2 6 9 D4 



if P e T^ 2 ' 1 ) (R 16 ), 2h = 0, 2, 4, 6, 8, 12, 
if Pe ft( 10 )(R 16 ), 
if Pe W( 12 )(M 16 ), 

if p = 1, 

ifPeH (2,i> (K 16 ), 2/i = 2,4,6, 

if p e w( 8 >(M 16 ), 
if Pe W( 10 )(M 16 ). 



Every shell of A is a spherical 7-design, and we obtain a cubature formula of 
strength 9 by combining two different shells of A, for example A4 and A@. We also 
obtain a cubature formulas of strength 11 of support -^(AuA')4U-^(AuA') 6 , and a 

cubature formulas of strength 13 of support ^(AuA')4U^(AuA') 6 U^=(AuA') 10 . 
6.4 The D 4 root lattice 

Let A be the lattice generated by the root system D4. It is a 2-modular lattice. 
Here, we will apply Lemma El with 

G = Aut(AuA'). 

We now describe A, A' and G. We have 

A = [x = {x\, ar 2 , X3, X4) | Xi € Z, X\ + x 2 + £3 + ^4 € 2Z}, 

A' = {x = (xi, x 2 , X3, X4) I Xi € ^75 Z, — £j € v2Z}. 

The group Aut(A) contains (actually, is) a reflection group, denoted by W(F4), 

■a, where a is one of the 



which is generated by the 24 reflections x 
following vectors: 



,{x\c_ 
' (a\a 



(±1,0,0,0), (±1,±1,0,0), (±1,±1,±1,±1) 

(all choices of signs and all permutations of coordinates). The invariant polynomi- 
als of this group are known Smit95, §7.4, pp. 217-218]; they form a polynomial 
C-algebra of basis {/i 2 , he, hs, /112}, where 

^2 = sym(x 2 ), he = sym^xfx^) — 3 sym(x 2 X2X 2 ), 

hg = sym(xf) + 14 sym^x^) + 168x 2 X2X 2 xf, 

hi2 = sym(xj 2 ) + 22sym(x^X2) + 165 sym^XjXg) 

+ 330 sym^x^XgX 2 ) + 330 syn^xfx^XgX 2 ), 
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and where 



sym(p) 



l(6 4 ) P | 



(64)2 := {er € ©4 | o ■ p = p}, 



cr£&4 



where the symmetric group 64 acts by permuting the coordinates. For our purpose, 
it is better to choose the basis 



Ho := h 2 



Hq := 8/i6 
The orthonormal transformation 



hi 



H u := 64/i 12 - 55h 8 hl ~ 176/ii + 220/i 6 /i| - ll/if. 



10&8 

'2 



7^, 



T 



/-V2/2 y/2/2 \ 

-a/2/2 

V2/2 V2/2 

V -V2/2 V2/2J 



exchanges A and A'. So, G = Aut(A U A') is generated by Aut(A) and T. The 
action of T on the polynomials Hk is: 



THo — Ho 



TH e , 



-Ha 



TH* 



H. 



THis 



-H12 



Therefore, the polynomials that are invariant by G are linearly generated by the 
polynomials 

H 2 H qH2H s 12 , (3 + 5 = mod 2. 

Let p(' £ )(R™) be the space of homogeneous polynomial functions on R™ of de- 
gree k. Since ft( fe )(R«) i s the kernel of the surjective G -equivariant Laplacian 
A : 7>( fc )(M") -> V {k - 2) (R n ), we have 

d G := dimft (fe) (R n ) G = dimV {k) (R") G - dimP^" 2 ) (R") G . 

We can compute 

1 + X 18 



(1 - X 8 )(l - A 12 )(l - X 24 ) 



= 1 + X + A + A ib + 0(X i8 ). 



So, we have by Proposition [7| 
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A,P 



e 



A',P 



'2 6 D4 ifP=l, 

if FeH^flf 

Cl (F)A 16 e D4 ifPeft( 8 )(R 8 ) G , 

^(PjAfae^ if P6HM^ G 



{0}, 2/1 = 2,4,6,10,14, 



Thus, using Lemma [§] together with Proposition we obtain cubature formulas of 
strengths 11 and 15 by combining two or three shells of A U A'. Note however that 
certain combinations of shells, such as -^(AU A') 2 U -^j(A U A') 4 are not possible, 
because Lemma[!5]can provide a degenarate linear equations system. (In the present 
example, we have indeed ^j(A U A') 2 = ^(Au A') 4 .) 

6.5 The shorter Leech lattice 

Let A be the unique (up to isometry) odd unimodular integral lattice of dimension 23 
and of minimum 3; it is the shorter Leech lattice. As for even ones, the theta series 
of odd unimodular lattices can be computed: see |Pach-a| . So, we can apply our 
method to construct cubature formulas from the shells of A. 
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7 Cubature formulas obtained by our construction 

The cubature formulas we obtain by our method are summarized in Tables lol to 1141 
These tables involved the lattices listed in Tables 01 and [5] 

In Tables ED to ^] the column "bound" indicates the lower bound of the size 
of cubature formulas of strength t on § n_1 obtained by Linear Programming. The 
indication in parentheses has the following meaning: 

• (LPd) refers to the estimate of the Linear Programming bound obtained by 
the procedure described in Section [3 paragraph "Numerical estimate of the 
Linear Programming bound" , where 2c? is the degree of the polynomial used for 
the estimate. The calculations have been performed with "Maple" software, 
using the "Optimization" package. Because of limitations in time and memory 
of our computer, we have not always chosen an optimal d, in which case we 
write (LP^g?); similarly, when the bound we have found can almost certainly 
be improved, we put the sign ^ before it. 

• (T) refers to the Delsarte or (Fisher-type) bound; it is the same as (LPd) 
when 2d + 1 = t (the strength). 

• When it is known that tight spherical designs with the corresponding param- 
eters (n, t) do not exist, we note (D) to indicate the Delsarte bound increased 
by 1 

The columns "set" and "shells" indicate respectively the set E and the shells 
uij used in the construction of Section according to the notation of Lemma 03 
The column "size" gives the size of the cubature formula so constructed. 

In the case n = 4, we have also listed in Table[7|the cubature formulas of smallest 
known sizes for fixed values of t; 

Most spherical designs with n ^ 24 described in BacVenOl, Table 3, p. 108] 
appear also in our tables. 



n 


£ 


notation 


name 


4 


2 




D4 root lattice 


8 


1 


r s 


Korkine-Zolotareff 


12 


3 




Coxeter-Todd 


14 


3 


Ql4 


Quebbemann 


16 


2 


BW 16 


Barnes- Wall 


20 


2 


N W jy(2) ^O) 
JV 20 > JV 20 ) JV 20 


Nebe 


24 


3 


N 24 


Nebe 


24 


1 


A 24 


Leech 



Table 4: Even ^-modular lattices involved in Tables IB1 to 1141 



n 


notation 


name 


23 


O23 


shorter Leech lattice 



Table 5: Odd unimodular lattice involved in Table 
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strength 


set 


shells 


size 


bound 


5 


D 4 


2 


24 


21 (D) 


7 


D 4 UD' 4 


2 


48 


42 (LP5) 


11 


D 4 \JD' 4 


2,6 


240 


120 (LP11) 


15 


D 4 UD' 4 


2,6,10 


528 


267 (LP22) 



Table 6: Cubature formulas for n = 4 — see also TabkEI 



strength 


best known cubature formula 


size 


bound 


3 


root system Af 


8 


8 (T) 


5 


root system D4 


24 


21 (D) 


7 


|HadSlo94j 


46 


42 (LP5) 


9 


announced in HadSlo94 


86 


74 (LP7) 


11 


vertices of the 600-cell 


120 


120 (LP 11) 


19 


vertices of the 600-cell and of 
its dual, the 120-cell ISali75l 


720 


500 (LP>18) 



Table 7: Best known cubature formulas for n = 4 



strength 


set 


shells 


size 


bound 


7 


r 8 


2 


240 


240 (T) 


11 


r 8 


2,4 


2 400 


1856 (LP8) 


13 


r s 


2,6,8 


24 240 


4 361 (LP12) 


15 


r s 


2,4,6,8 


26 400 


9190 (LP16) 



Table 8: Cubature formulas for n = 8 



strength 


set 


shells 


size 


bound 


5 


Ku 


4 


756 


157 (D) 


7 


K12 U K[ 2 


4 


1512 


729 (D) 


9 


K 12 U K[ 2 


6 


8 064 


2 940 (LP6) 


11 


K 12 U K' 12 


4,6,8 


50 400 


10 604 (LP7) 



Table 9: Cubature formulas for n = 12 
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strength 


set 


shells 


size 


bound 


5 


Ql4 


4 


756 


211 (D) 


7 


Qi4 u Q' 14 


4 


1512 


1121 (D) 


9 


Ql4 U Q'14 


4,8 


89 964 


4902 (LP6) 


11 


Qi4 u <2' 14 


4,6,8 


107436 


20 817 (LP7) 



Table 10: Cubaturc formulas for n = 14 



strength 


set 


shells 


size 


bound 


7 


BW 16 


4 


4 320 


1633 (D) 


9 


BW 16 


4,6 


65 560 


7 753 (D) 


11 


BWi 6 UBW' 16 


4,6 


131520 


37166 (LP7) 


13 


BW 16 UBW' 16 


4,6,10 


4 555 200 


^ 146153 (LP8) 



Table 11: Cubaturc formulas for n = 16 



strength 


set 


shells 


size 


bound 


5 


^20 


4 


3 960 


421 (D) 


7 


iV 20 U N' 20 


4 


7 920 


3 081 (D) 


9 


N 20 U 


4,6 


345 840 


17 711 (D) 


11 


N 20 U N 20 


4,6,8 


4 527 600 


95 309 (LP7) 



In this table, 7V 2 o denotes any of the three lattices N$ , N^q\ and Nl 
Table 12: Cubaturc formulas for n — 20 



strength 


set 


shells 


size 


bound 


7 


o 23 


3 


4 600 


4 600 (T) 


9 


o 23 


3,5 


958 458 


29 901 (D) 


11 


o 23 


3,4,6 


6 574 550 


166 808 (LP7) 



Table 13: Cubaturc formulas for n — 23 
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strength 


set 


shells 


size 


bound 


5 




6 


26 208 


601 (D) 


7 


N 2 a U 


6 


52 416 


5 201 (D) 


11 


A 24 


4 


196 560 


196 560 (T) 


15 


A 24 


4,6 


16 969 680 


^ 6179 991 (LP9) 


17 


A 24 


4,6,8 


415 003 680 


^ 27131261 (LP13) 


19 


A 24 


4,6,8,10 


5 044 384 800 


> 116 303 274 (LP15) 



Table 14: Cubaturc formulas for n — 24 
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